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We investigate several mechanisms of vortex creation during splitting of a spinor 
BEC in a magnetic trap controlled by a pair of current carrying wires and bias 
magnetic fields. Our study is motivated by a recent MIT experiment on splitting 
BECs with a similar trap, where unexpected fork-like structure appeared in the 
interference fringes corresponding to interference of two condensates, one with and 
the other without a singly quantized vortex. It is well-known that in a spin-1 BEC in 
a quadrupole trap a doubly quantized vortex is produced topologically by a "slow" 
reversal of bias magnetic field Bz- We find that in the magnetic trap considered 
it is also possible to produce a 4- and 1-quantized vortex in a spin-1 BEC. The 
latter is possible, for example, during the magnetic field switching-off process. We 
therefore provide a possible explanation for the unexpected interference patterns in 
the experiment. We also give an example of the creation of singly quantized vortices 
due to "fast" splitting, which is a possible alternative mechanism of the interference 
pattern. 
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I. INTRODUCTION 

Coherent manipulation of matter-waves is presently a very important experimental and 
theoretical field. Coherent splitting of matter waves into spatially separate atomic wave 
packets with a well-defined relative phase is necessary for atom interferometry and quantum 
information processing applications. Bose-Einstein condensates (BECs) open new unprece- 
dented perspectives in this direction. Coherent splitting of a BEC was recently realized both 
in an optical double- well trap |l| and in a magnetic chip-based double- well trap In the 
latter experiment, the magnetic trap was produced by a pair of current carrying wires and 
a bias field used to control the distance between the wells |2i]. An intriguing feature of this 
experiment was appearence of a fork-like structure in the absorption image of interference 
fringes designating creation of a singly quantized vortex in one of the wells. Thus perturba- 
tions during condensate manipulations were violent enough to generate vortices. To identify 
possible mechanisms of vortex creation is an inevitable challenge for future pursuit of BEC 
interferometry with such kind of an experimental configuration. 

Here we study two independent scenarios that could lead to vortex creation. One of them 
is the phase imprinting during the switching-off process (due to different decay times of the 
bias fields and the fields of the two wires), while the other is dynamical vortex creation 
during fast splitting. Although the switching-off time in the particular experiment ^] was 
only about 20 /is (a duration much shorter than the inverse of any trap frequency), it is 
slow as compared to the Larmor frequency in the magnetic field of order of IG. Therefore, 
adiabatic imprinting might take place. On the other hand, the authors of Ref. P| think the 
phase imprinting mechanism is an unlikely explanation for the vortex creation since they 
have never observed interference pattern of a doubly quantized vortex (which is created 
topologically in spin-1 condensates when zero point of the magnetic field crosses through 
condensates). We find realistic scenarios for a singly quantized vortex interference pattern to 
appear as a result of topological phase imprinting during the switching-off process. Indeed 
as we assume exponential decay of magnetic fields produced by the two wires and the 
bias field Bz with different decay constants, with parameters close to the experiment we 
obtain numerically that a large part of the atoms (about half of the condensate) can be 
transferred to the component \I^o (i-e-? with m = in z-quantized basis, see Eqs. (0) and 
(j2)) for notation) which has singly quantized vortices around each of the two minima of the 
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magnetic field. In contrast, when we try to obtain doubly quantized vortices using bias field 
Bz reversal, it turns out to be more difficult. At such short timescales, only a small part 
of the population is transferred to the component with doubly quantized vortices, while 
the rest of the population is redistributed among other components. Creation of doubly 
quantized vortices in such a configuration requires times of the order of hundreds of fis. We 
also present an example of dynamical vortex creation during fast splitting. In this case, 
the condensate (including the component) acquires winding phase during the splitting 
process. This leads to clear fork-like structures in all components and thus the total density 
upon expansion. Although in our calculations this requires a sufficiently shorter time of 
splitting than that which was actually effected in the real experiment Q] , it might be caused 
by the oversimplified model (without jitter or fluctuations of the magnetic field, etc). 

The rest of the paper is organized as follows. In Sec. II we introduce the model (spinor 
BEC described by 2D GP equation with spin degrees of freedom), describe the configuration 
of the magnetic field given by the two-wires setup, and its weak- and strong-field seeking 
states. In Sec. Ill, we present numerical results on vortex creation in the system. In Sec. IIIA 
we study several examples of how an unexpected topological phase imprinting could have 
taken place in the MIT experiment. There we take into account gravity and use realistic 
parameters relevant to the experiment. We assume that after switching off the magnetic field 
the field decays faster than transversal magnetic field Bj_, which can lead to the phase 
imprinting. In Sec. IIIB we present examples of dynamical vortex creation during splitting 
process. When split slowly in 30 ms, the condensate develops no vortex. However, reduction 
of the splitting duration to about 5 ms, vortices are created. Then, even in the case where 
the switching-off process leaves all population in the initial component, expansion of 
condensates produces very clear forks in the interference pattern. 



II. THE MODEL 



A. Spinor Bose condensate 

We consider a BEC of alkali atoms with hyperfine spin F = 1 using the z-quantized basis 
jj]. The order parameter is expanded as 

J2 ^m\m>, (1) 

m=±l,0 
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where |m > are eigenvectors of F^. 



1 >= (1, 0, 0)^, |0 >= (0, 1, 0)^ , I - 1 >= (0, 0, 1) 



(2) 



We use the time-dependent form of the GP equation with spin degrees of freedom developed 
in J 
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where (a = x, y, z) are the angular momentum operators in the basis of the eigenvectors 
of F.: 
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and Qn = ATi%^{ao + 2a2)/3m is the spin-independent and Qs = A'n%^{a2 — ao)/3m is the 
spin-dependent interaction coefficients. Here the scattering lengths ao and a2 characterize 
collisions between atoms with total spin and 2, m is the atomic mass. 

We consider a condensate of ^^Na, which was used in the experiment of j^. Parameters of 
Na are M = 3.81 • 10~^^ kg, a2 = 2.75nm, ag = 2.46nm 6|. Spin-independent potential ^(r) 
is provided either by gravity as in the MIT experiment {V{r) = —Mgy, we consider this 
case in Section 3), or, in principle, "optical plug" can be used. We assume the wavefunction 
of either component does not depend on the z coordinate, i.e. the condensate is quasi-2D. 
We therefore use two-dimensional nonlinearity parameters g"^^ and g^^ which are related as 
9n^/9s^ = 27.44. Configuration of the magnetic field is described in the next subsection. 



B. The magnetic trap 

We use a two-wires setup suggested in Ref. P]. It was used in a recent experiment 
of Y. Shin et al. jj], where several important elements were added making the system 
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essentially three-dimensional. We consider a two-dimensional trap here to simplify numerical 
calculations. Magnetic field produced by two parallel current-carrying wires is given by 



^ ^ ~^ d —d + x , . 

^ {x + dy + y^ {d-xY + y^' ^ ' 



where 2d is the distance between the wires. 

Additional magnetic field is added along z direction, and bias field B^ is added in 
X direction: B^ = B^ + B^ . For weak-field-seeking atoms, amplitude of magnetic field 
B{x,y) plays the role of trapping potential (while the potential for strong field-seeking 
atoms is —B{x,y)). At certain critical B^ = B^^ the trap potential B{x,y) has a single 
minimum located on distance d away from the surface, at the middle of the two wires. At 
B^ greater or smaller than Br^^ the potential has two wells: when AB^ = B^ — 5^.^ > 0, 
the two wells are separated in the x direction, while for AB^^ < they are separated in the 
y direction and has equal x-coordinates of their centers. Addition of bias magnetic field B^ 
along the y-axis rotates these two wells in the xy plane. It is important that the critical 
single well potential is harmonic (corresponding to hexapole configuration of magnetic field), 
while either of the separated double wells is quadrupole. We assume uniform B^ is applied 
to the system. In numerical calculations, we used the parameters close to the experiment 
of Ref. Q: d = 150/im (this corresponds to the distance between the wires of 300/im), 
= 24G, 5,(0) = IG. 



C. Strong- and weak-field seeking states in two- wires trap 

The important feature of static magnetic traps is that it confine the weak-field seeking 
state(s) (WFSS) which has a higher energy than the strong-field seeking state(s) (SPSS). 
Therefore, there is some difficulty in numerical preparation of initial states, since straight- 
forward imaginary time propagation would lead to SFSS. One need to search for a solution 
using a WFSS ansatz which is derived from the eigenvector of the B matrix. The eigenvalues 
of B are ±B ( where B = \J B^ + By + Bl) and (the latter corresponding to the neutral 
field-seeking state, NFSS), and the eigenvector corresponding to WFSS is 
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Suppose we have configuration with AB^ = (so the wells are coalesced in a single 
harmonic well). When [B^ ± iBy)/Bj_ = e^^"^ in the vicinity of the minimum of the well, 
where is the polar angle around the point of the minimum. The weak-field-seeking state 
order parameter has the form 
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where ip is the amplitude (common for all three components). At large B^, almost all 
population is in component. In order to avoid vorticity in this component, we should 
choose w = 2. Substituting ansatz (jHl) into GP equation, one obtains the one-component 
equation for the amplitude ip which can be solved using imaginary time propagation to find 
a ground state. However, it is also possible to propagate the original three-component GP 
equation in the imaginary time, restricting the solution to the anzatz (jH)). We used both 
algorithms in our calculations. 

One can see that one obtains a vortex with vorticity 4 from the initial non-vortex state in 
the critical hexapole trap by reversal of Bz- When AB^ 7^ 0, we have two quadrupole traps. 
In that case reversal of B^ leads to doubly quantized vortices in either trap, as was studied 
in several works recently (for example, Ref. j^). It is important that the component \l/o in 
that case contains singly quantized vortices. 



III. VORTEX CREATION DURING REALISTIC MAGNETIC TRAP 

MANIPULATIONS 



A. Phase imprinting 

Bias fields Bz and B^ and the two-wires magnetic field B^ are created by different 
sources. During a switching-off process, they can behave differently. Examples of nontrivial 
consequences of nonsynchronous decreasing processes of magnetic fields in other situations 
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were considered in Refs. 0, Ol- In order to achieve topological imprinting, we consider here 
the switching-off process with decaying faster than Bj_ (in our model calculations, we 
assume the fields and B^ decay synchronously, while the field B^ decays faster). We 
found a time of the order of tens of /is is enough to transfer a large part of the atoms to 
component \E'o with singly quantized vortices in each well. 

In contrast, we found that in order for B^ reversal to produce considerable population of 
\E'i component with doubly quantized vortices, much longer times are required, i.e. hundreds 
of fis. In the experiment of Ref. switching-off time was Toff ~ 20/is. Larmor frequency 
at -B = IG is ~ 700 kHz corresponding to tl = 1.4/is. Therefore, despite very short 
switching-off time, "adiabatic" phase imprinting still might take place. We checked this 
guess numerically. 

Four stages of the process were carefully examined: preparation of initial state, splitting, 
switching-off process, and the free expansion. Firstly, we prepare the initial state in the 
single (critical) harmonic well using imaginary time propagation, as described in the previous 
Section. (In the real experiment, the condensate was prepared initially at AB^ = — 140mG 
in the bottom well of the double-well potential, but it is not important for our present 
purpose). When, the condensate was split by ramping Ai?^ from to 100 mG in real time. 
During this stage, almost all the population is in \E'_i component, because the B^ field of 
IG is large as compared with magnetic field B± in the vicinity of the minima of magnetic 
field, where the condensate resides. Provided the splitting is slow, \I'_.i component remains 
without vortices (however, at fast splitting it acquires phase winding, this case is considered 
in the next subsection). During the third stage, the magnetic field was turned off. Faster 
decay of B^ may lead to creation of two singly quantized vortices in component \l/o in either 
of the two wells. This would result in two fork-like structures in the interference pattern of 
this component (only one of the forks points up, because charges of the topological vortices 
in either well are the same (positive)). The interference fringes were formed during the 
forth stage of numerical calculations: expansion of the two condensates without magnetic 
field. For numerical purposes, we also turned off the gravity field in the forth stage, but we 
suppose this does not infiuence results significantly. In the real experiment, typical time of 
expansion was about 22 ms. In our calculations, we typically calculate up to 4 ms only. 

In Fig.l, fast splitting of the condensate is shown. The condensate was prepared in the 
merged well (at AB^ = 0). Splitting was done approximately in 5 ms, much faster than in 
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FIG. 1: Splitting of the condensate by ramping AB^ from to 100 mG. Total density is shown. 
The origin of the coordinates is the merge point without gravity. Parameters are close to the 
experiment of Shin et. al; splitting was fast: time of splitting about 5ms. 

the experiment. In Fig. 2, slow splitting of the condensate is monitored (time of splitting 
about 30 ms). The final state is slightly different about 2 % of the initial population was 
lost because during very slow splitting. The atoms that make transitions from WFSS to 
SFSS fly away (they fly up towards the wires and are removed by absorbing potential on 
the boundaries of the mesh). The atoms that make transitions to NFSS, fly down because 
of the gravity (and are also removed). 

Densities and phases of the components of the split condensate are presented. Component 
has no vorticity, while \l/o has singly quantized vortices in both wells and has doubly 
quantized vortices. Since = IG, almost all population is in \E'_i component without 
vorticity. 

In Fig. 3, interference fringes formed during expansion of the condensates are shown. 
Magnetic fields were switched off in approximately 20/is in such a way that no transition 
from component occured {B^ was decreased more slowly than B±). 

Then, we consider a case where B^ is decreasing faster than In Fig. 4, the process 
of topological vortex formation is shown during switching-off process. Magnetic fields were 
decreasing exponentially, with Bz decaying faster, and were turned off completely at t = 
0.06r ^ 20/is (r ^ 350/is is the characteristic time period used in the program). About 
half of the total population were transferred to the component \E'o. An important feature 
of nonadiabatic transitions during process of switching-off can be noticed. The part of 
condensate in component residing far from the minima of the magnetic field is more 
easily converted to other components than the part near the minima. As a result, decreasing 
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FIG. 2: Slow splitting of the condensate by ramping AB^ from to 100 mG. Time of splitting 
about 30ms. For each component, density and phase profiles are depicted. In the component \I'o, 
singly quantized vortices of topological nature are seen. Their phase singularities reside in the 
minima of the magnetic field which slowly move as ABx is increased. 
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FIG. 3: Formation of fringes during expansion of the slowly split condensates. Switching-off of 
magnetic fields was done in such a way that almost all population (more than 99%) remains in 
^'-1 component. 
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FIG. 4: Formation of singly quantized vortices in component during switching-off process. -B^ 
decay faster than B±. Time of switching-off process is about 20 /xs. 

Bz slices out a part of the condensate around the minima of the magnetic field and leaves 
it in the initial \E'_i component. This is due to the fact that in the region where B_i is 
not zero, a nonvanishing gap between WFSS and NFSS remains; besides, the stronger the 
field B± is, the more slowly the total magnetic field rotates, making it easier for the spin 
of an atom to follow it. Quantitatively, nonadiabatic transitions due to nonsynchronously 
decreasing trapping magnetic fields (in other configurations) were considered recently in Ref. 
0] by generalizing the Landau-Zener formula for the multilevel case. It can be seen also 
that during fast switching-off process, the total density is almost unaffected, therefore the 
nonlinear interaction coefficient Qn plays no role in the process. 

Time evolution of populations of the components is shown in Fig. 5. In Fig. 6, interfer- 
ence fringes formed during expansion of the condensates are shown. Fork-like structure is 
seen in the component \l/o. 

Topological vortices in component \l/o in both wells has equal charges. Therefore, only 
one fork points up (Ref. ll|). The \l/o component rotates clockwise during expansion due to 
equal charge of the vortices, and the fork is moving up. Because of this, although during our 
calculations it has not appeared in total density profile, we believe that after a longer time 
it will show up (note that almost half of the total population is in the \l/o component). The 
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FIG. 5: Time evolution of populations of the components of the condensate during switching-off. 
Time of switching-off process T^fj is about 20 /xs (at this time, all remaining magnetic fields are 
abruptly turned off ). 

absence of the fork in the total density profile on early stages of expansion is related to the 
above-mentioned behavior of the switching-off process that leaves \E'_i component residing 
near the minima of magnetic field almost unaffected. After the switching-off this component 
therefore is concentrated at the same places where the phase singularities of component ^'o 
are located. 

For clarity, we give also an example of vortex formation without gravity (Fig. 9). Here, 
the both forks are clearly seen. 



B. Dynamical vortex creation 

In Fig. 7, fast splitting of the condensate is monitored. Dynamically created vortices 
appear in all the components (in component \E'o dynamical and topological vortices coexist). 
Centers of dynamically created vortices (the phase singularities) lie outside the condensate, 
but branch cuts go through the condensate causing phase winding across it. So this situation 
is different from the so-called ghost vortices appearing in studies of stirring condensates (Ref. 
j^). Although vortices are almost not visible in the density profile of the condensate, its 
expansion leads to characteristic fork-like structures in the interference fringes due to the 



12 



40 
20 

-20 



"(D « 

O Q) 
I- X) 









n 


n' 








IB 




CO 











-20 20 




20 

X (|i m) 



(b) 



5.000E-4 



(a) 



FIG. 6: Formation of fringes during expansion of the condensate with imprinted vortices. Almost 
half of the population were transferred to the component ^'o during switching-off. Topological 
vortices in component ^'o has equal charges in both wells. Only one fork points up in the \I'o 
component. In total density pattern, the fork have not appeared yet. However, note that in 
the experiment expansion took much longer time. We see that the \I'o component with vortices 
rotates clockwise during expansion. After a long time, the fork therefore might appear in the total 
density pattern too (our current computational recourses forbid such a long-time propagation), 
(a) Density of components ^'o,^'-i, and total density, (b) Density of the component with the 
fork- like structure after approximately 4 ms of expansion. 
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FIG. 7: Fast splitting of the condensate. Dynamical formation of vortices: vortices are formed 
in all components (note that more than 99% of total population is in component ). In the 
^'o component dynamical and topological vortices coexist. Topological vortices are the same as in 
Fig.2. 

phase winding. 

Magnetic field was switched off in such a way that almost no transfer to \l/o occur. So, \l/o 
component does not influence the dynamics. In Fig. 8, interference fringes resulted from the 
subsequent expansion are shown. Two fork-like structures (which appear in all components) 
are clearly seen. Vortices in the two condensates have opposite charges, therefore the forks 
point up (splitting of the condensate on two parts in the presence of gravity effectively causes 
rotation of each part in opposite directions). 

To finish off, let us notice that the calculations reported were done using the sin-DVR 
method (with a spatial mesh of 200 x 160 grid points), and time propagation using the 
Runge-Kutta method (the method is the generalization of the one used in Ref. [l^ to the 
case of spinor condensates). We try to model the system as close to the experimental one 
as possible. However, we found that, for example, Qs factor does not affect the results, i.e., 
the results are qualitatively the same with = 0. The most important deviation of the 
model from the experimental system is its reduced dimensionality. The real system is three- 
dimensional, with nonuniform magnetic fields, with an additional two pairs of wires creating 
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FIG. 8: Formation of fringes during expansion of the rapidly split condensates. During switching- 
off, almost no transfer occur to the component ^q. Dynamical vortices in the two wells have 
different charges, therefore the two forks point up. 

magnetic fields to compensate partly for nonuniformity of the base magnetic fields and 
asymmetry, etc. Such complication is presently beyond our available resources. However, 
we believe the study of the idealized model do allow to discuss different mechanisms of 
vortex creation in the experiment. 

IV. CONCLUSION 

The authors of the experimental work P| supposed the phase imprinting mechanism to 
be unlikely for explaining the appearance of the fork-like structure and that "the observed 
phase singularity definitely shows the breakdown of adiabaticity" . To the contrary, we found 
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FIG. 9: Dynamics without gravity, (a) Formation of singly quantized vortices in ^'o component 
during switching-off process. Bz decay faster than B±. Time of switching-off process was about 
20 /is. (b) Fork-hke structures in the density of ^'o component. 

a realistic scenario based on nonsynchronous decreasing processes of the magnetic fields can 
explain the phase singularity even within the assumption of adiabatic evolution. 

However, we found that fast splitting can also lead to dynamical vortex creation, and 
that the dynamically created vortices produce interference patterns with the forks of better 
contrast. Detailed study of these processes is left for future research. In any event, the 
splitting of a spinor BEC is a rather violent process so that further consideration of atom 
interferometry with spinor BEC is necessary. 
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